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Black branes are studied in Einstein-Gauss-Bonnet (EGB) gravity. Evaporation drives black 
branes towards one of two singularities depending on the sign of a, the Gauss-Bonnet coupling. For 
positive a and sufficiently large ratio ^/a/L, where L/2tt is the radius of compactification, black 
branes avoid the Gregory-Laflamme (GL) instability before reaching a critical state. No black branes 
with the radius of horizon smaller than the critical value can exist. Approaching the critical state 
branes have a nonzero Hawking temperature. For negative a all black branes encounter the GL 
instability. No black branes may exist outside of the interval of the critical values, < /3 < 3, where 
P = 1 — 8a/r^ and rn is the radius of horizon of the black brane. The first order phase transition 
line of GL transitions ends in a second order phase transition point at /3 = 0. 

PACS numbers: 04.50.+h, 04.70.Bw, 04.70.Dy 

INTRODUCTION 

Gravity in higher dimensions has been in the forefront of research for a considerable time. Most of the studies start 
with the conventional Einstein action generalized to D > 4 dimensions. In the low energy limit quantum theories 
of gravity like string theory yield additional higher order curvature correction terms to the Einstein action. [T| Some 
higher order curvature terms when taken in isolation and quantized can potentially lead to problematic high energy 
behavior, such as the presence of negative norm states. The simplest term avoiding these problems is the Gauss- 
Bonnet (GB) term. In this letter we analyze the critical behavior of a class of classical solutions to Einstein gravity, 
augmented by the GB term. Such a term arises in the gravity action in string theory when corrections to zero slope 
limit are calculated [5] [3]. For simplicity, we restrict our discussion to five space-time dimension, but as we will point 
out later the results can be simply generalized to higher dimensions. 

In Einstein gravity D — A solutions, like a black hole can be trivially extended to a black string {D — 5) or black 
brane solutions {D > 5). Solutions become more complicated if higher derivative terms are added to the Einstein 
action, such as the GB term. Though asymptotically Minkowski black hole solutions have been found in D > 5 
theories |3] , no exact black string (brane) solutions are known in EGB gravity. Yet the investigations of black strings 
(branes) is more important because with compact extra dimensions all but the lightest static objects must be black 
strings (branes). 

Two alternative techniques have been used to investigate black brane solutions. Kobayashi and Tanaka [4^ solved the 
5 dimensional Einstein equations modified by the Lanczos tensor (variation of the Gauss-Bonnet term) , numerically. 
Using an expansion around the event horizon they also found an exact lower bound for the black string mass, Mc ~ y/a, 
where a is the coupling constant of the GB term. For AI < Mc the solutions do not have a horizon, they represent 
naked singularities. 

In a subsequent work [5^ , following |4| , we used a series expansion in a to investigate black strings in D > 5 EGB 
theory. We found that in every dimension, D > 5, a lower bound, similar to that in D = 5, exists for the mass of 
the black string. We also investigated the thermodynamics of the black string but the fifth order expansion employed 
in our paper was insufficient to get a definitive answer for the thermal behavior of the system when the mass of the 
black string approached the lower limit. 

In this paper we also use improved numerical and expansion methods to investigate black string solutions including 
the neighborhood of critical points, Pcrit — and 3, where the dimensionless GB coupling is defined as 

/3=l-8^. (1) 

' h 

In ([ij r/i denotes the radius of the horizon. (3 — determines the critical radius. The critical radius has a one-to-one 
relationship with the critical mass, which is determined by the condition (3 > 0. 

In Section 2 we will discuss the uniform black string, which should have a higher entropy than a black hole, when 
the mass is large enough. In fact, for sufficiently high mass it is the only static state of the system. We will use a 
horizon expansion to elucidate the nature of two critical points /3 = and 3. We will concentrate on Z? = 5 but, 
except for numerical details D > 6 theories are not different. Then we will use numerical calculations to investigate 
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the properties of black string at large distances from the horizon. We discuss their thermodynamics and calculate 
their entropy. 

In Section 3 we will investigate the extension of the Gregory-Laflamme instability [6] and the notion of the non- 
uniform black string, which is the preferred state if the radius of horizon goes under a critical multiple of the radius 
of compactification [7j [8j [9j. We will also plot the phase diagram of black strings in the (3 — rh plane. 

Section 4. summarizes our results. We will also discuss the possible directions for the continuation of this research. 



UNIFORM BLACK STRING 



In this section we will discuss the extension of a 4-dimensional solution to D. If we assume that the metric gab 
is an asymptotically flat solution of the 4 dimensional Einstein equation then the trivial extension of this metric, 
obtained by defining gai = and g^ — 1 satisfies the D dimensional Einstein equation. Here, and in what follows, 
D > i > 5. We tacitly assume that the coordinates, labeled by i are compactified on circles, in which space we will 
maintain a (S^)^^* symmetry throughout this paper. However, the trivially extended metric is not the solution of 
Einstein- Gauss-Bonnet equation 

Gab = Rab — ^dabR + OiLab — 0, (2) 

where a is the coupling, the Lanczos tensor is 

Lab = —-^gab {R^ — ^R^'^Rab + if ) + 2RRab — ^RacRb ^ '^RacdeRlf'^ + '^RacdbR'^'^ (3) 

and K is the Kretschmann scalar, 

K = RabcdR'""''. (4) 

In an appropriate coordinate system the metric components of uniform black brane solutions depend only on a 
radial coordinate. Nonuniform black branes and the Gregory-Laflamme instability will be discussed later. The metric 
of nonuniform black branes is a periodic function of the compactified coordinate, w, as well. 



Horizon expansion and scaling near critical points 

Throughout this paper we will use dimensionless coordinates, scaled by the radius of the horizon, r/j, of the uniform 
black string. In particular, we choose the horizon variable, p defined by 

P= : (5) 

to replace the radial coordinate, r. When the GB coupling, a, is also rescaled, by the introduction of the parameter, 
(3 (See eq. [T]), then all reference to the radius of the horizon is eliminated from the equations of motion. The only 
parameter remaining in the equations of motion is (3. True values of physical quantities are restored by multiplying 
by a power of r^, required by their dimensions. 

Defining the metric of a uniform string, made dimensionless by factoring out rf^, as 

ds^ = -f{p)df + ^dp^ + {p+ 1)2^^(2 + h{p)dw^^ (6) 

f(p) 



we can solve the equations of motion in a power series of p (horizon expansion) . The series coefficients depend only 
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on p. Displaying three leading orders only, the expansions have the following form: 

r( s _ r ( ~^ + ^+^^";^ 2 -55/33 - 76/3^/2 _^ 61/32- 416/33/2 + 135/3 -8V;g+ 147 ^ 
"'^^^ ~ 4(3 - /3) ^ 72(/3 - 3)2/33/2 



+ 1 7 - 8/3 + /3' 
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2 4(3-/3)73'^ 
(1 + ^/P) (-147 + 231x//3 - 261/3 + 161^3/2 ^ 47^2 _ 35^5/2 + ^3 _^ 3^7/2) 



2 I 

72(/3- 3)2/33/2 ^ 
, (, 4(V?-1) ^ 2 (/32 + 4/33/2 - 16/3 + 47?+ 7) , , ^ 

The two constants, /i and /iq, are not determined by the horizon expansion alone. They depend on /3 only. They 
must be introduced in /(p) and h(p) to assure that the metric components are asymptotically Minkowski. They are 
determined in the process of numerical integration of the equations of motion. 

The expansion coefficients have two singular points: (3 = and (3 = 3. Note that for a < /3 > 1. The physical 
domain is < /3 < 3. The expansion coefficients are complex for /3 < and there are singularities outside the horizon 
for /3 > 3. Questions have been raised whether negative values of a are admissible [13]. Nevertheless, for completeness 
we will also investigate that region, as well. Near the singular points the convergence radius of the expansion shrinks 
to zero. The leading singular terms of the expansion coefficients (subscript k implies the coefficient of p'') are the 
following: 



fk ^ 


. ^3/2-fe^ 


(3 




gk ^ 


. /31/2-fc, 


(3 




hk 


^ ^3/2-fc 


(3- 


-(3)-K 



(8) 

Since the expansion coefficients become complex at /3 < the solution does not exist for negative /3, contradicting 
the assumption of the existence of a horizon or/and a regular expansion around the horizon [4j. This result can, 
however, be the result of a poor choice of the coordinates. Therefore, we also calculate scalars, like the Ricci scalar 
(i?) and the Kretschmann scalar {K), at the horizon. Being scalars, they must be independent of the scaling factors 
/i or /lo which can be gauged away by rescaling the coordinates t or w. Indeed, we obtain 

^ l + V?-3/3 + /33/2 
(3-/3)(l + V3) ' 
69 + 2^//3 - 13/3 - 12/33/2 _ ^2 _^ 2/3^/2 + /33 

(3 - /3)2(1 + V?)2 ^> 

We can draw several conclusions from ([9]). First of all, R and K are non-analytic at /3 = and 3, ruling out mere 
coordinate singularities. Furthermore, both of these constants become complex for /3 < 0, therefore no uniform black 
string solution with regular horizon can exist in that region. /3 < implies < V^cx and, as we will see later, r^, 
though not quite proportional, but is a regular and bounded function of the ADM mass. Thus, the (3 > bound 
imposes a lower limit on the mass of the black string. 

Furthermore, while all curvature invariants are finite at /3 = they are infinite at /3 = 3. This points to the 
possibility of the existence of a critical solution at (3 = but not a,t (3 = 3. We will discuss the /3 = critical solution 
in the following section. 



The critical string 

To investigate the critical solution we start with a rearrangement of the horizon expansion series. Keeping only the 
leading singular terms in every order of p we are lead to functions of a scaling variable y = p/ (3. Thus, one can define 
alternative expansions in series of scaling functions 4>i{y) , "/i{y) , and rjiiij). The scaling expansion is defined as 

f{p) = p |i + ^['^i(y) + VpUv) + ••■]} , (10) 
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with similar expressions for g{p) and h{p) in terms of series of scaling functions, ^i{y) and r]i{y). Substituting (10 1 
into ([2]) we obtain differential equations for the scaling functions defined in ( |To| ). The equation for 0i is 

[62/3(2y0i - 1) + 6/0'i]0'/ + 27y\4'[f + 6y\l2y^i - + 2(182/20? - 12y0i - 7y) = 0, (11) 

while 7i and rji can be simply expressed by as 71 = 201 + y(j)[ and 771 — 40i/3. Though we could not find the 



solution of ( 1 1 1 in analytic form the general behavior of 0i can be analyzed fairly easily. Its power series in y is 
provided by the most singular terms of horizon expansion ([t]) . Analyzing an 80th order series expansion we were able 
to ascertain that 0i has a singularity inside the horizon of the form 

o\3/2 

2/+ 7) • (12) 



Eq. (12 1 suggests that the metric becomes singular inside the horizon, at p = —3/3/7 and cannot be continued to 
smaller values of p. Note that this singularity moves to p = at /? = 0. We will indeed see below that the critical 
solution is singular at p = and cannot be continued to p < 0. There is no impediment, however to continue 0i to 



positive values. In fact the analysis of Eq. (Ill shows that the asymptotic behavior of 0i is 



^,^±yiy-'f' + Oiy-'). (13) 

Numerical integration of Eq. (Ill shows that for a function with the series expansion implied by ([7| the positive sign 
should be chosen in (13). Now in the (3-^0 limit the first terms of (10 1 become 

P 



/(P) = . 1. ^ 




P l + (14) 




The /? — > limit of higher order expansion terms lead to higher powers of p, including all half-integer powers. This 
can be ascertained by repeating the horizon expansion ([?]) with half-integer powers included (starting from p'^^^). This 
leads to the following expansion of the critical black string {(3 — 0) solution: 

J^-p) _ ^ I „ I ^./' „3/2 „2 ^^oao ^5/2 , 250322 ^3 

9ip) 

^ V o 

Kp) = ho (^l-^+^^l^p'^-' + ^p' + ...]. (15) 

If one attempts an expansion with half-integer powers of p for general (3 then either all the half-integer powers vanish 
or /3 = 0. Conversely, no horizon expansion with integer powers exists for the critical black string. 

The value of scalars R and K can also be calculated at p = using ([T5|. One obtains R — —2/3 and K = 92/3, 
precisely the values obtained from the (3 = limit of ([9|. In fact, one can show that all scalars are finite in the critical 
state. Correction to scalars, when one moves away from the horizon are of 0{y/p). Further evidence for the fact that 
the critical string is the limit of non-critical strings when /3 ^ comes from numerical integration of the equations of 
motion, to be discussed later. 

The p = surface cannot really be called a horizon, just the limit of horizons when P 0. While it is an infinite 
red shift null surface, the metric cannot be continued inside the p = surface. We believe that there is no way to 
facilitate such a continuation. Strictly speaking, though this is a string solution, it is not "black." It is however a 
limiting solution of a series of black string solutions in every quantity that we investigated. At all /3 > the horizon is 
completely regular and there is always a region inside the horizon, p > —3/3/7 , where the solution can be continued 
without any impediment. 

Time-like geodesies can be continued to the horizon both in the (3 > and /3 = cases. A test particle reaches the 
horizon in finite proper time, r. However, while in the /3 > case the test particle passes through the horizon into a 
region from which it cannot escape any more, for f3 = geodesies have a pathological behavior. Choosing a geodesic 
for a particle falling towards the horizon, p(t) < 0, and the time scale such that p(0) = the radial coordinate turns 
complex for r > 0. Therefore, the role of the critical solution as a valid state of black strings is questionable. This 
question will be discussed later in more details, following the discussion of numerical integration of the equations of 
motion and of the Gregory-Laflamme instability. 
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Numerical integration of the equations of motion 

We build on the results of the previous sections to obtain the metric functions at arbitrary p > 0. Before starting 
the numerical integration from the horizon we calculate the value of the metric functions and their derivatives from 
the horizon expansion at p = A > 0. This step is necessary because p = is a singular point of the differential 
equations. A is chosen to be small enough so that the last terms of the 30th order horizon expansion is smaller than 
10^^''. Then we integrate the equations numerically from p = A using metric ([6|. The step size, S, in the numerical 
integration is chosen such that S < O.OIA. 

Expansions ([7| contain undetermined constants /i and Hq. In an initial run these constants are chosen to be 1. 
Then, though the calculated metric functions tend to a finite value as p oo, their asymptotic limit is not 1. Fitting 
a form 

+ ^ + ^ + (16) 

to /(p), and a similar form to h{p) we define /i — l/af and ho — 1/ah to get, after a second round of numerical 
integrations, metric functions that are asymptotically 1, which corresponds to Minkowski space. 



Thermodynamics 

The numerical form of the metric functions allows us to investigate the thermodynamic properties of the black 
string. The first law of thermodynamics states that 

dM KdL 

dS^— —, (17) 

where M is the Arnowit-Deser-Misner (ADM) mass, T is the Hawking temperature of the black string, k is the 
tension, and L is circumference of the compact string. 

The surface gravity on the horizon, k, is related to the Hawking temperature as follows [10] [11] [12] 



Using (|6| and ([t]) we obtain 



k=^^, (19) 
r/. 2^5(0) 



where we restored the correct dimension of the surface gravity by reinstating a factor of r^^ ^ . As we pointed out earlier 



/'(O) and 5(0) are determined in the process of numerical integration. They are functions of (3 only. Combining (18 1 



and ( 19 1 we can write the Hawking temperature as 



T.— W). (20) 

Form factor Ft(/3), normalized to -Ft(I) = 1, is calculated in the process of our numerical integration procedure. It 
is a smooth function of (3 between /3 = and /? = 3. It tends to a finite value at /3 = 0. That limit, along with the 
limits of the other two form factors defined below, agree with the values obtained from the critical string at (3 — 0. 
Ft {(3) — > cx) in the limit /3 ^ 3 implying that T — s- 00 as /? ^ 3. Ft and the other two form factors are plotted in 
Fig.l, 

The ADM mass and the tension are expressed simply by the long range components of /(r) and h{r) [9 . Namely, af- 
ter renormalizing /(p) and h{p) so that /(oo) = 1, h{oo) = 1 and defining the leading, 0(p^^), asymptotic coefficients 
as /_i and h-i, then 

r = ^(/-.-2/._.) = ^^.(/?), (22) 
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where Fm {P) and Fr (/?) are two additional form factors depending on /3 only. They are both bounded on the interval 
< /3 < 3. These form factors also converge smoothly to the corresponding values for the critical string at (3 = 0. 
They are plotted in Fig 1. along with Ft{I3). All three of the form factors are normalized to 1 at /? = 1 (vanishing 
GB coupling). The three form factors are not independent from each other. The intcgrability condition 




FIG. 1: The plot of the scaling functions Ft{[3) (solid line), Fm{P) (dotted line) and -FV(/3) (dashed line) as a function of /3 
between the two critical points, and 3 



dL 



dM 



relates them. (23 1 implies 



(23) 



(24) 



(1 - f3){2FM - Fr)Ffr + Ft {-Fm + + (1 - (3)f;) - 0. 

The calculated values of the form factors satisfy (24l to 0(10^"*), except near /3 = 3, where the rapid variation of Ft 
makes the evaluation of the numerical derivative difficult. 
We can also calculate the entropy. We have 



S- 



dM 



87ra 



where q{L) is independent of M. (25) can be rewritten by changing variables as 



8£7ra f Fm + 2(1 - (3)F 'j, 
G J {I-P^Ft 



M 



G 



liL) 



(25) 



(26) 



The function q{L), a yet undetermined arbitrary function of L, can be fixed if we take the derivative of (26 1 with 
respect to L and use the intcgrability condition ( 24 1 . Setting the lower bound of the integration at /? = we obtain 
the final form of the entropy for a > as 



SLna 
G 



Fm(/3) + 2(1 - P)F'M ^ ^Fm{0) - FriO) ' 



(1 - /3)2Ft(/3) 



Ft(0) 



(27) 



An overall constant was fixed in (27 1 so that 5 = for L = 0. From our numerical integration we obtain [2Fm(0) — 
JV(0)]/Ft(0) = 2.05614. 

For completeness we write down an equivalent expression for a < 0. Considering that Ft diverges around the 
critical point /3 = 3 the integral over f3 is convergent at that point so we can write 



S=- 



SLna 
G 



fM(/9)+2(l-/3)n,(/3) 
(1 - I3YFt{P) 



dp, 



(28) 
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implying that the entropy vanishes at the critical point, /9 = 3. 



The limit of S when a can be calculated if we change the integration variable in (27 1 to r/j. Then we obtain 
2L7r 



S 



G 



Fm{PK)) + r'^dFM{PK))ldri^,^ ^ ^J^^'^^^ " ^^^^^ 



Ft{Q) 



Since I3{r'j^) = 1 — ^ajr'^ and the form factors are equal to unity at /3 = 1, at that point (29 1 reduces to 



5*0 



2L7r 



Lnrl 



(29) 



(30) 



which is the expression for the entropy of a black string in Einstein gravity. The entropy can be represented in the 




0.5 1.0 1.5 2.0 2.5 

FIG. 2: Plot of function s(/3) 



3.0 



fori 



(31) 



at all < /3 < 3. The function s{j3) obtained from numerical integration is plotted in Fig. 2. Note that s(l) = 1 and 
s(3) = 0. 

GREGORY-LAFLAMME INSTABILITY AND NONUNIFORM BLACK STRINGS 
A simple picture of the life of a black string 

We will investigate what happens to black strings as Hawking radiation reduces their mass. Consider the ADM 



mass as a function of Vh and j3, (21 1. A/ is a monotonic function of r/j, because for the whole physical range of (3 



dM 



2G; 



.{Fm + 2{1-(3)F'm)>0. 



(32) 



With decreasing rh (3 decreases if a > and increases if a < 0. Typical trajectories of black strings in the (/?, 2Trrh/L) 
plane are shown in Fig. 3. According to ([IJ if a > > y/Sa. When hits the value \/8a the singular point, 
/3 = 0, is reached. In the critical state the black string has a finite Hawking temperature, but it cannot decay and 
stay in the black string state any more. We will return to the discussion of its further fate later. 

If a < /? increases with decreasing r^. Eventually, it will hit the critical value P — 3, where — \/—Aa. However, 
the Hawking temperature diverges at /3 = 3, so approaching this point the black string will completely evaporate. 
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Before turning to the investigation of instability we briefly discuss the fate of black branes in D > 5 EGB theories. 
We will use the notation d = a/rf^. The coefficients of the horizon expansion depend on d only. Their singularity 
structure in /3 is the same as that at D = 5, except /? has a more complicated expression 

f3 has a single zero for a > 0, which zero varies with D from a minimum of a = 0.1233 at Z? = 6 to a = 0.1443 at 
D = oo. For negative a we have /3 > though /3 is not a monotonic function of d. Still, just like for D — 5 the 
horizon expansion coefficients are singular at a = —1/4. Thus, aside for small numerical changes the fate of black 
branes is the same as that of black strings. They are driven to the critical state /3 = if a > and to a = —1/4 if 
a < 0. 



Gregory-Laflamme instability 

In the absence of the GB interaction term (a — 0), uniform black strings become unstable if 27rr/i/L < k ~ 0.876 [5]. 
At this point a nonuniform mode becomes marginally tachyonic, a zero mode. As it was shown by Gubser [8J, Harmark 
and Obers, Kol, Piran, and Sorkin ^ this point corresponds to a first order transition to a stable nonuniform black 
string state. Such a state is associated with a periodic dependence of the metric function on the compactified 
coordinates, in our case w. It is fair to assume that such an instability also affects the EGB theory when a ^ 0. 
Therefore, the simplistic picture about the life of a black string which we depicted in the previous section may be 
misleading. The black string may hit an instability line before it reaches either of the singular points (/3 = or 3). 
Therefore, we performed a stability analysis of EGB black strings using our numerical approximation to the metric 
components of the uniform string solution. 

Following Refs. jH] and [S] we perturbed the metric with a periodic function with wave number k of of the compact 
coordinate, w. Considering the rescaling of coordinates with k can be identified with n = 2'Krh/L. Applying the 
Landau-Ginzburg theory of phase transitions to the Einstein black string it was found fS] [SI that at the Gregory- 
Laflamme point [B], K ~ 0.876, the system undergoes a first order phase transition to a nonuniform state with a 
finite non-uniformity. It is quite reasonable to expect that a first order transition exists at nonzero values of the GB 
coupling, as well. 

We consider the following static perturbation of metric Q [S] [5] 

ds^ = -e^^f{p)df + e^^^dp^ + {p+ Ife^^'dn^ + e^"h{p)dw^, (33) 

J[P) 



where the exponents A, B, C, and H arc functions of both p and w. Though in (33 1 we gauged away a possible 
non-diagonal term still the gauge has not been fixed completely. A perturbation analysis is applied to ( 33 1 , assuming 
a small deviation from the uniform string. It is also assumed that the smallest possible wavelength, k dominates the 
perturbation. Higher order corrections, along with higher harmonics can be considered in a systematic perturbation 
analysis that we leave to future work. Accordingly, we use the following ansatz for the perturbed metric 

A = 1 + X a{p) cos{kw) . (34) 

We also introduce amplitudes 6(p), c{p) and h(p) for the perturbative form of coefficients B, C, and H. Linearization 
of the equations of motion in A leads to a set of equation for the amplitudes a, b, c and h. These equations still have a 
substantial gauge freedom. As we prefer to work with a single amplitude, we use the Kol-Sorkin gauge [9J to eliminate 
b, c, and h, using the gauge freedom. 

Making the final gauge choice c{p) — —a{p)/2, Grw — can be solved for b{p). Here and in what follows we use 
the notation G^i, for the components of the equation of motion. In the next step Grr = and Gee — are used to 
eliminate h{p) and its derivatives. Then Gu = and G^w = provide identical second order differential equations of 
the form 

W{p)a" + X{p)a' + [Y{p) + K^Z{p)]a = (35) 

for a(p). 

The coefficients W , X, Y, and Z are extremely complicated functions of the metric components /, g, h and their 
derivatives |14j . An important feature of W is its proportionality to /. Consequently W vanishes on the horizon 
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making (35) singular at p = 0. The ratios X/W ^ p , Z/W ~ 1 and Y/W ~ p ^ when p — > oo. Thus, (35 1 is a 



Schrodinger-like equation, with —k as a bound state eigenvalue. Our aim is to determine the single eigenvalue of the 



equation for every choice of (3, which is the only parameter entering (35) 



Due to the singularity at the horizon we cannot integrate (35) starting from the horizon. However, we generated a 



30th order Taylor series of a around p — 0, using the known series expansion coefficients of the metric components of 
the uniform black string (/, g, h). This allows us to calculate a and a' at p = A, which is the same quantity which 
has also been used at the numerical integration of the equations of motion for the uniform black string. 

We use the shooting method, starting from p = A to find the eigenvalue. At each f3 the coefficients of (35) are 



calculated using the stored values of functions /, g, and h and their derivatives. We determine k at a large number of 
values of f3 between the two singular points, and 3. The solid line of Fig. 3. represents the eigenvalue k = 27rr^/L 
as a function of (3. While at /3 = k reaches a finite value, Kc, it diverges at /3 = 3. This implies that for a < black 
branes must undergo a Gregory-Lafiamme transition before reaching the critical value and evaporate completely. The 
limit of K as calculated from /3 7^ solutions agrees with k calculated from the critical solution. 




FIG. 3: Phase diagram of black strings. The solid line represents the first order transition line at 2nrh/L = «:(/3) as a function 
of /3 between the two critical points, and 3. The dotted, dash-dotted, and dashed lines represent time-lines of black strings 
at different values of a 



The time-lines of black strings are depicted by dash-dotted, dotted, and dashed lines in Fig. 3. These time-lines 
depend on a only. At fixed a different black strings may start their life at different points but they follow the same 
line. If 

r{a, L) = > Kc = 2.2868 (36) 

then black strings completely avoid the instability transition before reaching the critical state (dash-dotted line in 
Fig. 2). However, if r{a,L) < Kc then every black string undergoes an instability transition and turns into a non- 
uniform black string before it could reach the critical line. The trajectory of black strings at such a value of a is 
represented by the dotted line. Note that r{a, L) is completely determined by the geometry of space and the GB 
coupling. Furthermore, if a < then /3 > 1 and the time-line approaches /3 = 3 (dashed line). Finally, if a = then 
black strings follow the thin vertical dashed line at /3 = 1. 



Non-uniform black strings near /3 = 

In contrast to uniform black strings non-uniform black strings require the solution of a two variable problem. 
Therefore, it is quite a bit more complicated to investigate the fate of non-unifom black strings. However, we can 
substantially simplify matters if we restrict ourselves to sufficiently small values of /3. Though we cannot quite follow 
non-uniform black strings through the region of instability, we assume that if a > they eventually also approach the 
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/3 = line in the plot of Fig. 2. Of course, it is quite possible that nonuniform black strings also undergo an instability 
transition, before they could reach the critical state. Still the following calculations, pertaining to the behavior of 
nonuniform stings near (3 = 0, certainly apply to a neighborhood of the critical point, /? = 0, 2Trrh/L ~ Kc- 



The simplest way to proceed is to investigate the horizon expansion (33 1 with u;-dependent coefficients. We write 
the following ansatz for the metric 

ds' = -(p+"^p' + ..) dt'+^-^^^'^'l^!^^^^^^ (37) 
\ 2 ) p+ -fp^ + ... 

Note that the w dependence can be gauged away from gu and from the zeroth order contribution to gww Expand- 
ing the equations of motion in p, the lowest order contribution to Grw implies that the periodic function 6o(if) is 
independent of w. In the next step one can solve the lowest order contributions to Ggg = and Gww = equations 
for hi{'w) and hi{w). Then the lowest order contributions to equations G^r = and 0^,^ = are identical. They 
are second order equations for ci{w), containing cq{w) and its second derivative. Solving the resulting equations for 
ci(w) the discriminant of the equation must be positive for a real solution. This discriminant has the form 

d = 16[co(u;)]4 + (1 - mc'oH]' - 8(1 - /3)[co(«^)]'[2 + 4{w)] - (1 - /3)2co(u;)[c^(H]'[2 + i4{w)] 

+ 2(1 - f3)[co{w)r Hc',{w)r + (1 - f3)c'^{w)[2 + ^{w)]} . (38) 

Now we know that for uniform black strings co{w) = 1 and for a small non-uniformity co{w) must have the form 

co(w) — 1 + Acos(ku'). (39) 

Then expanding c? in A and f3 and keeping the leading order contributions only, we obtain 

d~ 414:13+ {4 + K^)Xcos{nw)]. (40) 



( 40 1 implies that the discriminant is positive for all w only if 



4+ K'^ 



First of all (411 implies that A — > when we approach the critical state, P = 0. This makes our approximation of 
expanding d in A correct for small /?. Secondly, this implies that the non-uniform black string becomes a uniform 
black string again at the critical line. Thirdly, this implies the line of first order Gregory-Laflamme transition to 
the non- uniform state ends in a second order transition point at /3 = 0. This is a well-known feature of first order 
transitions. 

CONCLUSIONS 

We have studied uniform and non-uniform black strings in Z) = 5 EGB gravity. We paid particular attention to 
the singularities at the cndpoints of range, i?, < /3 < 3, where (3 is defined in ([T|). No black string solutions exists 
outside R. In particular, we found a critical solution at /3 = which, in contrast to solutions at /? > 0, has a horizon 
expansion that includes half-integer powers of the radial horizon variable, p ([5]). Not only the metric components, but 
all scalars, including the Ricci scalar and the Kretschmann scalar, are also singular (though finite) at /3 = ([9|. 

At P > the metric components are singular inside the horizon, at p = —3/3/7 ( |12[ ). It is always possible to continue 
the metric to a region inside the horizon. Introducing Eddington-Fonkelstein coordinates 



the metric takes the form 

rfs^ = -f{p)du^ + 2^/^du dp+{p+ ifdn^ + h{p)dw 



2 



regular everywhere including a neighborhood of the horizon. |15| The singularity at — 3/3/y moves to the horizon when 
/3 — » 0. We have not been able to construct a metric inside this singular point. At /3 = 0, the critical metric cannot be 
continued inside the horizon. This can be shown if near the horizon we transform the metric to Kruskal coordinates. 
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We define the Kruskal coordinates for the critical string such that the coefficient of dU dV is —1. They can be obtained 
in a asymptotic expansion from the Eddington-Finkelstein coordinates u and 

v = t+ fdp^^ 



fip) ■ 
We obtain 

ds'^ = -~dU dV + {p+ ifdn'^ + h[p]dw'^ 
where p has the following expression by the Kruskal coordinates U and V 



p^ ~uv{i + cVuv + ...), 

where c is a non- vanishing constant. It is obvious we cannot cross the light like surfaces [/ = or y = which 
represent the horizon. We also showed that a geodesic of a massive probe falling toward the horizon reaches it in 
finite proper time (at r = 0), but it cannot pass the horizon because its radial coordinate turns complex at t > 0. 
Therefore, the critical solution does not represent a black string. 

We used numerical techniques to calculate the solutions for the whole range R. We calculated the ADM mass, 
tension. Hawking temperature, and entropy. These quantities are represented by form factors, depending on P only, 
measuring the deviation from the a = (/9 = 1), Einstein gravity solutions. These form factors are defined in 



( 20 1 ( 21 1 ( 22 ) ( 31 ) .They are plotted in Figs 1 . and 2. While the Hawking-temperature tends to a finite value at /9 = 0, 
it diverges at /3 = 3. The ADM mass and tension are modified only by a finite amount compared to pure Einstein 
gravity and bounded from above and below over the whole range, R. We gave a closed expression for the entropy in 
terms of an integral over the form factors, which was evaluated. The entropy is a monotonically decreasing function 
of /3, vanishing in the limit (3 — 3. This is not surprising in view of the divergence of the Hawking temperature. 

We paid particular attention to the life line of black strings decaying via Hawking radiation. We found that every 
black string follows a trajectory in the (/?, k = 2-Krii/ L) plane, approaching (depending on the sign of a) one of 
the singular lines in the (/3 — r/j)-plane. This is depicted in Fig. 3. However, some of the black strings, depending 
whether 2it^/%cx/ L > Kc = 2.2868, will encounter a line of first order singularities, which is the extension of the famed 
Gregory-Laflamme singularity. This line is the solid line shown in the phase diagram of Fig. 3. Under this line only 
non-uniform black string are stable. As we commented before, non-uniform black strings may become unstable at a 
subsequent phase transition line. 

Finally, we showed, by investigating non-uniform string solutions near the horizon that the amplitude of non- 
uniformity vanishes near the /3 = line. This implies that non-uniform strings become uniform once again as they 
reach the P — line. As a corollary, we concluded that the line of first order phase transitions ends in a second order 
transition point at (3 = 0. 

There are two unresolved problems concerning the fate of black strings. If a > they are driven toward the critical 
state at /3 = 0. That state, however is very problematic. Though gtt has a linear zero and all scalars are finite at that 
zero, test particles that reach the p = point in finite time have nowhere to go. because the p(r) becomes complex 
as a function of proper time. 

As finding solutions for black holes in compactified spaces is even more difficult than finding black string solutions 
we have no way of comparing the entropy and relative stability of these solutions. This is one of the reasons why we 
are not able to answer what happens with black strings that reach the end of their lives at the (3 = line. Since the 
limit of their Hawking temperature is 

T, = Ft(0)-^^0.85--^, (42) 
47rV8Q; 4-77 V 8a 

not zero, as a study based on l/p expansions lead us to believe earlier, so they are certainly not frozen at that point. 



Their ADM mass is (21) 



M. = Fm(0)^.1.23^. (43) 

They cannot go into a black string state with a smaller ADM, mass, however, because such black string states do not 
exist. It seems that the only possibility is that some other, yet undiscovered, type of higher entropy static or possibly 
dynamic state exists, into which the black string can morph. Another possible solution, at least for £) > 6, is that 
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the inclusion of higher order Lovelock terms would cure this anomaly. The fate of black stings at and beyond /3 = 
is one of the most important questions we would like to investigate in the future. Other questions include a possible 
numerical study of black holes in compactified spaces. 

This work used a Kaluza-Klein type compactification. It would be more difficult, but perhaps more interesting to 
repeat this work in a Randall- Sundrum type ADS space. 
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